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We discuss simple generic model of “jet quenching” in which matter absorption is defined by one 
parameter. We show that as absorption grows, the azimuthal asymmetry parameter V 2 grows as 
well, reaching the finite limit uj which has a simple geometric interpretation. We show that this 
limit is still below the experimental values for 6 > pt > 2 GeV, according to preliminary data from 
STAR experiment at RHIC. We thus conclude that “jet quenching” models alone cannot account 
for the observed phenomenon, and speculate about alternative scenarios. 


1.Azimuthal asymmetry for non-central heavy ion col¬ 
lisions have been predicted [Q to be larger at RHIC then 
at lower energies. In hydrodynamic models this happens 
due to the stronger push by high pressure of Quark-Gluon 
Plasma well above the phase transition region, which is 
expected to be produced at RHIC. In contrast to that, 
models based on string picture of hadron production (e.g. 
RQMD and UrQMD event generators) or on mini-jet sce¬ 
narios (e.g. HIJING) have predicted its decrease. The is¬ 
sue has been settled already by the first data from RHIC, 
by STAR collaboration , which have found large asym¬ 
metry consistent with hydrodynamic predictions. De¬ 
tailed studies have provided significant details, such 
as the asymmetry parameter 

V 2 =< cos{2(j)) > (1) 

(where </> is the angle between the impact parameter 
and momentum of a secondary hadron in the transverse 
plane) as a function of centrality, particle type and its pt ■ 
Data from STAR and PHENIX experiments for rather 
wide range of momenta pt < 2 GeV agree well with these 
predictions for all secondaries. This is probably the most 
direct signature of QGP plasma formation, observed at 
RHIC. 

In this letter we however discuss a different ques¬ 
tion, related with (much less certain experimentally) such 
asymmetry at higher transverse momenta 6 > pt > 
2 GeV. According to the latest STAR data [|| (which are 
still unpublished and are thus considered preliminary, al¬ 
though reported at many meetings), V 2 {pt) for all charged 
secondaries seem to be about constant, for each central¬ 
ity. This means a different regime seem to be established 
in this region of pt, and the original intention of this 
note was to compare these data with a simple geomet¬ 
ric model for jet quenching by relating the asymmetry to 
the strength of the jet quenching itself. However, after 
playing with different versions of the model, from more 
complex to a most generic one to be reported in this note, 
I concluded that the intended fit is simply impossible. 

A “jet quenching” idea has been discussed for a long 
time, see e.g. [||, and it has been naturally related to 
the azimuthal asymmetry for non-central collisions. If 
a high-pt jet is loosing energy in matter, jet emission is 
dominated by the surface of the almond and the correla¬ 


tion between position and the emission direction appears, 
thus the observed azimuthal asymmetry. 

A relation between this phenomenon and data has 
been discussed in |^], where it was concluded that a 
combination of jet quenching and hydrodynamical ex¬ 
pansion can approximately describe them^ Later STAR 
data have shown at high p( = 2 — 6 an approximately pt- 
independent V 2 , which disagree with a decreasing trend 
expected from jet quenching. Qualitative discussion of 
many possible scenarios which can have such a behavior 
has been made in Ref. §, including the interplay of jet 
quenching, hydrodynamical expansion and “baryon junc¬ 
tion dynamics”. We return to this discussion at the end 
of the paper. 

2.The present work ignores such details as pt depen¬ 
dence of V 2 and focuses instead on its measured values: 
we demonstrate that looking at pure geometric aspect of 
the problem one can show that those are too high for any 
jet quenching model (without hydro). 

The most generic model we use can be described as 
follows. First, the distribution of origination points for 
outgoing jets is simulated: this is done using the usual 
assumption of parton model and the simplest model of 
nuclei as two homogeneous colliding spheres. (Diffuse 
boundary only makes effects smaller.) 

The second step is the calculation of the chances for 
the parton to escape the absorption in matter, as it goes 
out of the almond. The absorption rate is characterized 
by one (and the only) free parameter of the model. Its 
magnitude determines the strength of jet quenching it¬ 
self (the fraction of escaping partons f{pt,b)), with the 
predicted azimuthal asymmetry, V 2 {b). 

As high-pt partons move with the speed of light, we 
ignore possible change of shape due to geometrical ex¬ 
pansion of the “almond” during this time. (If anything, 
this will reduced the asymmetry, as expansion reduces 
spatial asymmetry.) 

Naturally, in the absence of an absorption there is no 


* Although no curves for jet quenching alone are shown, the 
text implies that it is indeed insufficient by itself, in agreement 
with the (more general) argument we will give below. 
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azimuthal asymmetry, U 2 (k = 0) = 0, while increasing 
absorption creates increasing V 2 - Interestingly, in the 
limit of very strong absorption the asymmetry reaches 
a finite limit, denoted by asterix below 

V2{b, K ^ Qo) V2{b) ( 2 ) 

The reason for that is that in this case all the emitted 
partons/hadrons originate from the thin surface of the 
almond (see below). Even in this case, however, partons 
have half solid angle open for them: thus (6) has direct 
geometric interpretation. The main point of this letter 
is that, after evaluating V 2 ib) values for the experimental 
conditions and comparing it with data we have found 
that even the limiting ones, v^ib), are below the data. 



FIG. 1. A frontal view of two colliding nuclei, with defini¬ 
tion of the axis. The black dot (x,y) inside the almond is the 
origin of the parton, which propagates in the direction of the 
unit vector n. 

3.Let us now provide more details about the model 
itself. In fig.l we show geometry of the collision and 
definition of two longitudinal lengths, L±(x, y) for a hard 
collision at point (x,y). For hard spheres 

L±{.x,y) = 2[R^- {x±b/2YY/'^ (3) 

The probability of production of a parton in hard colli¬ 
sions at position x,y is simply proportional to the prod¬ 
uct of longitudinal lengths P(x,y) = aL+{x,y)L-(x,y). 
Vanishing of each of these factors defines the boundary 
of the initial almond in the transverse plane. Fig.l shows 
the sketch of the initial distribution in transverse, x- 
y, plane. One characterizes it by the standard spatial 
anisotropy 

S2{b) =< - x^ > j <y^ P x^ > (4) 

where angular brackets means average over all produced 
jets, with the weight given by the parton model as de¬ 
scribed above. The distribution depends on impact pa¬ 
rameter b, indicated in the l.h.s. In the table below we 
will make integration over b with geometric weight 2TTbdb 
over bins of centrality, within limits defined by upper and 
lower percentage of the total cross section. 

The probability to escape depends not only on the 
point of jet origin but also on the optical depth of matter 


along the outgoing line, {x + s * Ux, y + s * Uy) which we 
calculate as follows 

nOO 

f = exp[—K / ds{L-L+){x + s * nx,y + s * Uy)] (5) 
Jo 

The parameter k (dimension fm~^) includes both the 
density of the material and the absorption rate. The 
following fig.2 shows how the efficiency of the parton 
quenching and the V 2 parameter depend on it, in the 
whole dynamical range. The dependence of jet quenching 
and V 2 on the absorption strength is shown in Fig. 2(b). 
It displays the saturation of V 2 as well as the tendency to¬ 
ward the surface emission at large absorption, mentioned 
in the introduction. 

The main outcome of the simulations is summarized 
in the Table 1, in which we compare the high absorp¬ 
tion limit V 2 calculated from the model with STAR pre¬ 
liminary dataQ Q at pt > 3.5 GeV. As one can easily 
see, even in the high absorption limit the model fails to 
reproduce data, being systematically below the present 
preliminary data. The difference is especially striking for 
the most central bin, in which the observed V 2 nearly 
matches the asymmetry S 2 of the original almond. 

(Although the result is described as that in the high ab¬ 
sorption limit, it actually corresponds to the calculation 
in which absorption was large but finite, k = .2fm~^, 
with the actual quenching factors / also given in the ta¬ 
ble.) 


Centrality % 

</> 

Vils2 

* 

V2 


0-11 

.018 

.32 

.042 

.12± 0.02 

11-34 

.027 

.35 

0.12 

.16± 0.02 

34-85 

.046 

.31 

0.16 

.22± 0.02 


TABLE I. The limiting momentum/spatial asymmetry for 
three different centrality selections of STAR, given as V 2 ver¬ 
sus the percentage of total AuAu cross section. The quantity 
< / > is the escape probability (5) averaged over produced 
jets in the collisions, with all directions and origin points. 


^The error bars are calculated by the author, based on three 
STAR points at the largest pt bins, for each centrality. As at 
this point the data still have preliminary status, the reader 
should be warned that the error bars may be modified and the 
systematic errors be better understood and included. Now it 
is not possible to quantify the problem we discuss any better. 
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FIG. 2. Dependence of (a) the strength of “jet quenching”, 
the fraction of escaping partons /(ft) and (b) the correspond¬ 
ing asymmetry parameter V 2 (k) on the absorption strength 
parameter k, in fm~^ . Figure (c) show the (un-normalized) 
distribution over initial position x of the escaped partons, 
for the cut \y\ <2 fm. The solid and dashed lines are for 
ft = 0.2, 0.05 fm~^ respectively. All three figures correspond 
to the most peripheral centrality bin, 34-81 percent. 

5.Let us now summarize the main result of this let¬ 


ter: the dynamical range of “jet quenching” scenarios is 
approximately confined in the region 

0 < v*2/s2 < 1/3 (6) 

while the preliminary STAR data give larger values 
V 2 IS 2 = .5 — 1 and therefore they cannot be explain by 
models of this kind alone , no matter what magnitude of 
jet quenching be used. The generic model used can of 
course be modified in many ways, but it seems unlikely 
that jet quenching by matter absorption in whatever form 
is able to explain these data by itself. 

Assuming these preliminary STAR data are correct, let 
us consider what their explanation can be. The main 
shortcoming of the model comes from the idea that sec¬ 
ondaries in this region of pt originate only from jets, 
obtaining azimuthal asymmetry only from geometrical 
asymmetry of the almond. The interplay between jet 
quenching and hydro expansion, quantitatively discussed 
in Q , only reduces the effect due to reduction of geomet¬ 
rical asymmetry with time. 

The resolution of this puzzle can only be obtained if a 
signihcant fraction of secondaries originate from a source 
other than jets. A general discussion in Q have men¬ 
tioned a possibility that V 2 for baryons and pions can 
be very different, with the former getting a contribution 
from “baryon junctions” and/or collective flow, as the 
sources complementary to jets. We also think that it is 
likely to be the explanation, although we are quite scep¬ 
tical about the role of the baryon junctions^]. 

Collective hydro expansion is not just a simple and gen¬ 
eral concept, it is basically the only known mechanism 
capable to generate very large values of the azimuthal 
anisotropy. (Let us remind the reader why is it so. Due 
to different hydro motion in different directions, spectra 
have the (/-dependent pt slopes, resulting in asymme¬ 
try, V 2 , about linearly increasing with pt .) However, 
the issue is far from being simple, and a significant role 
of hydro component in the high-p( tails of spectra, at 
Pt ^ 4 — 6 GeV, is a very non-trivial thing. These tails 
of the particle spectra are 6 orders of magnitude below 
the majority of the particles, way below where a macro¬ 
scopic language is routinely used. More work is needed 
in order to understand whether such approaches can at 
all be used in this region. In connection with that let me 
mention a very interesting paper by Molnar and Gyu- 
lassy ^ in which V 2 has been generated kinetically in 
some model with very large cross section, way above per¬ 
turbative predictions. Although it is far from being clear 
that the extreme assumptions made in these calculations 


* Recent STAR data on spectra of </ mesons have provided 
one more argument against it. These data show that (p has 
Pt slopes consistent with hydro predictions ^ with the slope 
not very different from the nucleon’s: so it is the mass not 
the baryon number which matters here. 
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are realistic, it has been able to yield collective flow and 
sufficient values of the V 2 - 

Experimentally it is quite obvious what one should do: 
as soon as statistics will allow, to study at such pt for 
any identified secondaries. Particles which are seen via 
decays, e.g. A and Kg and especially (j) can be identified 
at rather high momenta and are thus most interesting. 
Since jets decay into pions much more than into strange 
mesons like (j) and especially into baryons, one should 
expect the corresponding fractions of jet-originated and 
hydro-originated secondaries be very different for all of 
them. The observed constancy of V 2 with pt for all charge 
secondaries is likely to be just a result of occasional can¬ 
cellation between rising hydro-based and decreasing jet- 
based components. 

Note Added in proofs After the paper was submit¬ 
ted to PRC, STAR data used have passed necessary pro¬ 
cedures and are no longer preliminary. The final data 
are submitted to PRL for publication, as “Azimuthal 
anisotropy and correlations in the hard scatt ering regime 
at RHIC”, by C.Adler et al, nucl-ex/020600^ . Systematic 
effects due to two-body correlations have been studied by 
comparison between 2-particle and 4-particle cumulants. 
When the latter values for V 2 is used, the discrepancy 
with the maximal model values at strong jet quenching 
V 2 nearly disappears. Another significant fact reported 
in this STAR publication is the first direct observation 
of jet component in the 2-body correlations. More recent 
STAR data at 200 GeV/N have been presented at re¬ 
cent Quark Matter 2002 conference. Due to much higher 
statistics, those extends to larger pt = Q — 12 GeV, but 
display about the same 02 - Good agreement between the 
measured values of V 2 and the theoretical high quench¬ 
ing limit V 2 has been shown in the summary talk there 
by S.A.Voloshin: it seem to suggest that geometrical in¬ 
terpretation of azimuthal asymmetry suggested in this 
paper seem to be correct, after all. 

This work is partially supported by the US-DOE grants 
DE-FG02-88ER40388 and DE-FG03-97ER4014. 
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